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Abstract 

Event-by-event analysis of heavy-ion collision events is an important tool for the study of 
the QCD phase boundary and formation of a quark-gluon plasma. A universal feature of phase 
boundaries is the appearance of increased fluctuations of conserved measures as manifested 
by excess measure variance compared to a reference. In this paper I consider a particular 
aspect of EbyE analysis emphasizing global- variables variance comparisons and the central limit 
theorem. I find that the central limit theorem is, in a broader interpretation, a statement about 
the scale invariance of total variance for a measure distribution, which in turn relates to the 
scale-dependent symmetry properties of the distribution.. I further generalize this concept to 
the relationship between the scale dependence of a covariance matrix for all conserved measures 
defined on a dynamical system and a matrix of covariance integrals defined on two-point measure 
spaces, which points the way to a detailed description of the symmetry dynamics of a complex 
measure system. Finally, I relate this generalized description to several recently proposed or 
completed event- by- event analyses. 

1 Introduction 



With the advent of high-multiplicity heavy-ion collisions at the CERN SPS and RHIC there is the 
possibility to extract statistically significant dynamical information from individual collision events. 
By determining the degree of correlation of the multiparticle final state for individual events within 
a large ensemble of events one can probe the dynamical history of collisions and the QCD phase 
boundary. This process is called event-by-event (EbyE) analysis. 

A central issue for EbyE analysis is the question whether or not collision events with similar 
initial conditions experience the same dynamical history, excepting finite-number fluctuations and 
known hadronic variance sources. If similarly prepared events differ significantly from one another 
beyond expectation, then analysis must proceed to categorize events into dynamical classes or on 
a continuum, and to provide a physics interpretation for nontrivial variations. 

Based on very general arguments and the observed behavior of normal bulk matter near phase 
boundaries we expect to observe significant changes in fluctuations near the QCD phase boundary. 
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The amplitude and character of these changes and their observabihty given finite coUision systems 
and non-equihbrium colhsion trajectories are major questions for EbyE physics. Experimentally, 
the problem calls for extracting all available information from each event by a complete topological 
analysis in order to obtain maximum sensitivity to dynamical fluctuations. 

Various measures can be used to categorize events by their 'morphology.' The common element 
of all measures is the correlation content of the multiparticle final state. Correlation measures 
can be loosely grouped into global measures, some of which have counterparts in thermodynamics, 
and scale-local measures, which explicitly determine the scale dependence of correlations. Given 
an event-wise measure we also require a corresponding comparison measure to enable quantitative 
event comparisons. 

In this paper I emphasize EbyE analysis techniques which relate to variances and second mo- 
ments. These global statistical measures are directly related to two-point correlation measures 
defined in various pair spaces. Higher moments are related to general g-point correlations and are 
more efficiently dealt with by a generalized system of scale-local EbyE analysis which will be treated 
in subsequent papers. The emphasis on variances and two-point correlations in this paper is both 
timely and provides a reasonably accessible introductory treatment of EbyE analysis. 

The paper describes the general types of EbyE analysis, the available system of correlation 
measures, the system of ensembles and its relationship to topological partition systems, considers 
maximally symmetric systems and the statistical model as archetypal correlation reference systems, 
the relationship between fluctuations and correlations, introduces the central limit theorem as an 
archetype for variance comparison, considers a fluctuation comparison measure (^pj as a difference 
between rms fluctuation amplitudes, analyzes the structure of as a specific variance comparison 
measure, introduces the concept of scale-dependent covariance matrices as a generalized measure 
system, relates scale-dependent covariance to a system of two-point correlation integrals, and con- 
siders the novel connections between scale invariance and the central limit theorem which arise 
from this work. 

2 Correlation Measures and EbyE Analysis 

Correlation measures provide the basic analytical infrastructure for EbyE analysis. A complete 
characterization of the correlation content of the multiparticle final state should exhaust all signif- 
icant information in the data. Comparison of event-wise correlations in a data system with those 
of a reference system (which may invoke 'equilibrium' or an assumption of maximum symmetry) 
offers a precise means to detect deviations from a model or ensemble average. Significant deviations 
provide a basis for discovering new physics. One then attempts to understand the physical origin 
of these deviations, leading to a revised model as part of an oscillatory (but hopefully converging) 
scientific process. 

2.1 EbyE analysis types 

We can distinguish two types of event-by-event analysis as represented in Fig. ||. In analysis type A 
each event is fully characterized for information content and recorded in an 'event spectrum,' a space 
designed to represent differing event morphologies in an efficient manner. This could be as simple 
as a 1-D space recording the event-wise slope parameter of mt distributions or the event- wise <pt>. 
Or, it could be a complex scale-dependent representation of correlations in (77, (p) distributions for 
a minijet analysis. Through the intermediary of event-spectrum distributions object and reference 
event spectra are compared for nonstatistical differences in event morphology. In this analysis we 



EVENT-BY-EVENT ANALYSIS AND THE CENTRAL LIMIT THEOREM 



3 



are usually not interested in the 'DC aspects of correlation common to every event, but only in 
those 'AC aspects that change from one event to another. 

Type A analysis generally requires large acceptance and large event multiplicity. One attempts 
to characterize all correlation aspects of each collision event over most of phase space, at least 
implicitly using g-particle correlations, where q may range up to a significant fraction of the event 
multiplicity. Questions of interest include large-scale variations in collision dynamics (dynamical 
fluctuations) such as flow fluctuations, baryon stopping fluctuations, jets, departures from chemical 
and thermal equilibration and chiral symmetry variations. The identity of individual events is 
preserved for essentially all aspects of this analysis. 



In analysis type B elements of each event are accumulated in an inclusive sibling space {e.g., 
a pion pair space for HBT analysis or an azimuthal space with registered event planes for flow 
analysis) and compared to a mixed-element space (no siblings by construction, hence minimally 
correlated) to determine differential inclusive symmetry properties of the event ensemble. Here we 
are interested in the 'DC aspects of correlation which are true for almost every event, but which 
may deviate from a model expectation. The identity of individual events is typically discarded in 
type-i? analysis. 

If statistical measures (such as moments) are extracted from an event spectrum obtained in a 
type-A analysis in such a way that event identities are lost then the result is effectively a type- 
B analysis, even though the statistics were obtained through a type-yl analysis. For instance, 
the moments of a < pt > distribution form the basis for a type-B analysis, although the < pt > 
distribution itself is the event spectrum for a type- A analysis. The variance of a <pt> distribution 
contains information which is closely related to that of an inclusive pair-correlation analysis. As 
we shall see however, this is a complex relationship. In contrast, an analysis in which features of 
an event spectrum are used to select special event classes is definitely a type-yl analysis. Special 
event classes can then be used as the basis for conventional inclusive analyses {e.g., particle ratios 
or baryon stopping) which would not be feasible on a single-event basis due to limited statistical 
power but may be essential to interpret the special event class. 

For type-_B analysis large acceptance and large event multiplicity are still highly desirable. For 
N-N vs A-A comparisons type-i? analysis may be the only practical study method due to the low 
multiplicity of N-N events. Comprehensive analysis of two-point densities over a large acceptance 
requires simultaneous measurement over the entire two-point space to minimize systematic error. As 
we shall see, in-depth EbyE analysis is often statistics limited and depends on access to the largest 
possible multiparticle phase space, practical conditions that require the largest possible detector 
acceptance. EbyE analysis based on variance comparisons as discussed in this paper represents 
only the most elementary initial approach. More complex extensions to arbitrary g-tuple analysis 
and scale-local analysis of higher-dimensional measure distributions will extend the initial studies 




correlation analysis 



Figure 1: Signal processing analog 
of EbyE analysis showing AC {A) 
and DC {B) aspects. Type A anal- 
ysis emphasizes the differences in 
correlation content among individual 
events. Type B analysis emphasizes 
the correlation content common to 
an ensemble of events relative to a 
reference population. 
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of simple 1-D distributions described here. 



2.2 Correlation measures 



The basic tools of EbyE analysis are correlation measures. Each multiparticle event distribution has 
a certain information content relative to a reference. In principle, all significant event information 
can be extracted to an efficient system of correlation measures by a data-compression process. If a 
distribution nears maximum symmetry within constraints the event information may be represented 
by a very simple measure system. We can divide correlation measures roughly into global measures, 
which are effectively integrals over some scale interval, and scale-dependent or scale-local measures. 
A central aspect of EbyE analysis is comparison between object distributions and references. Thus, 
every correlation measure requires an associated comparison measure. A summary of correlation 
measures is presented in Table |l]. 



measure 


comparison measure 


global measures 


moments 


central-limit tests 


Wm{D) =<x">i/"^ 




model distributions 


GoF criteria 


fix, a) 


x'{fix,a),D} 


scale-dependent measures 


factorial moments 


normed factorial moments 


H^g){D,6x) 


FqiDi,D2,6x) 


SCA 


Sq{D,6x) entropy 


Iq{Di, D2,Sx) information 


dq{D,6x) dimension 


Adq{Di, D2,Sx) dim. transp. 


Vq{D,6x) volume 


Pg {Di , D2 , 6x) density 



Table 1: Correlation mea- 
sures and associated compar- 
ison measures. D is a distri- 
bution or data set, {Di,D2) 
is an object /reference pair of 
distributions or data sets, a 
is a set of model parameters, 
GoF is goodness-of-fit, g is a 
g-tuple index, 6x is a scale 
(bin size) parameter. 



Global measures include distribution moments and model parameters. The respective compari- 
son measures are moment differences (e.g., central-limit tests) and goodness-of-fit (GoF) measures 
(x^ or likelyhood). The reference may be represented either by a set of moments or by a restricted 
region in the parameter space of a model distribution. Event- wise information then consists of 
deviations from the reference represented by moment differences or parameter deviations. As we 
shall see, global measures are effectively based on scale integrals, although this may not be obvious 
from typical usage. Integral measures are more appropriate in the case of low statistical power (low 
event multiplicity) and small correlation content, when only one or a few significant parameters 
can be extracted from an event. 

Scale- dependent measures become more important with increasing event multiplicity and/or 
correlation amplitude. It is then possible for an event to carry more information, and one can 
determine the scale and space dependence of correlations. A system of measures based on the 
scale-dependent correlation integral has been developed in two closely-related manifestations. Fac- 
torial moment analysis has a substantial history in multiparticle analysis Q and is based on the 
factorial moment as scale-dependent measure and the normalized factorial moment or moment ratio 
as the comparison measure. Scaled correlation analysis (SCA) is based on a system of generalized 
entropies [Q, ^ extended to scale-local form. The three closely-related topological measures entropy, 
dimension and volume are indexed by g-tuple number and scale. The corresponding comparison 
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measures are information, dimension transport and density. In some applications model distribu- 
tions such as the negative-binomial distribution have been employed in scale-dependent analyses 
IQ, ^. It is possible to make a connection between scale-dependent model parameters and SCA 
measures. 

3 Ensembles 

Many-body dynamics, equilibration and fluctuations - issues central to EbyE analysis - can be 
studied with an ensemble approach. We encounter ensemble and partition concepts repeatedly 
in this paper. An ensemble in the sense of statistical mechanics is a set of nominally identical 
dynamical systems (events) prepared with the same constraints but otherwise uncorrelated. Under 
Gibbs' ergodic hypothesis the phase-space distribution of system states of the ensemble at any time 
should reflect the phase-space history of an individual system state over a long interval. For nearly 
symmetric systems this hypothesis provides an essential reference. 

A generalized ensemble can be formed from an arbitrarily large bounded volume containing 
a measure system in 'statistical equilibrium' (maximally symmetric within global constraints). Q 
is invested with a set of conserved global extensive quantities (additive measures) such as total 
energy £, total volume V and total object number M. This list extends to all quantites conserved 
by a given dynamical system and is represented below by m = {rrii}. The global volume is 
partitioned arbitrarily into elementary volumes Ve- Each volume element 14 contains local samples 
of extensive quantities such as energy Eg and number A^e- Each sample set {E,N)e provides 
local estimators {E /N ,E /V,N /V)^ for global intensive quantities (T,e,p) representing the global 
constraints {£,V,Af). The (infinite) set of elements from all possible partitions of I designate 
the general ensemble (GE). This description of a general ensemble adheres closely to the language 
of Borel measure theory so that topological analysis methods can more easily be connected to the 
statistical study of dynamical systems. A system of ensembles is illustrated in Fig. ^. 

The measure contents of GE partition elements can be represented in a single-point measure 
space Gi{m) (cf Sees. ^ and |l^). The measure space in Fig. § is a subspace Gi{E,N). For a 
minimally correlated system the GE measure distribution is localized near the main diagonal of 
Gi(m). The degree of localization (fiuctuation amplitude or variance) depends on the effective 
total number of dynamical objects M in the space O and the fractional volume or scale of typical 
ensemble elements measured by position along the main diagonal. Note that fluctuation amplitudes 
decrease (are 'suppressed') as this partition scale approaches the total volume V of the space fi. The 
distribution becomes sensitive to the global constraints or boundedness of the dynamical system. 

The classical microcanonical, canonical and grand canonical ensembles are conventional subsets 
of the GE. The grand canonical ensemble (GCE) consists of all partition elements with some fixed 
volume V (system open to particle and energy exchange). The canonical ensemble (CE) is a special 
case of the GCE in which object number N is held fixed (system open to energy exchange). The 
microcanonical ensemble is a further special case in which {E, N, V) are all held fixed (closed statis- 
tical system), itself can be seen as an element of a microcanonical ensemble for a larger system. 
More generally, given an arbitrary measure system defined on a partitioned multidimensional space, 
one can define a hierarchy of sub ensembles of successively lower dimensionality corresponding to 
increasing degrees of constraint on partition elements and measures. 

For the purpose of EbyE analysis of nuclear collisions the GCE for an equilibrated system can 
serve as a minimally correlated (maximally symmetric) reference for the transverse phase space 
of an event ensemble. One asks whether the stopping process (source of total Pt) for a fixed 
macroscopic collision geometry is the same for each event, excepting finite- number fluctuations. 
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Figure 2: Ensembles: An isolated volume O, with extensive properties {£,N ■ ■ ■) is defined. Various 
subsets of the class of all possible partitions of can be identified as ensembles. Contents of 
partition elements are represented as points in the single-point measure space G\{fn). Points 
representing the general ensemble are localized near the main diagonal of this space. The projection 
of Gi{fh) for the GE onto the Gi{E,N) plane is shown as the shaded region in the figure. The 
grand canonical ensemble is represented by a slice through this distribution at fixed volume V. The 
result is a localized distribution (concentric circles) centered at {E, N) with variances and a]^. 

and whether the 'stopped' energy is fully equilibrated among transverse degrees of freedom. 

To answer these questions one must compare observed measure distributions for an event en- 
semble to a suitable reference such as the GCE for an equilibrated system. Multiplicity can be 
compared to a Poisson, quantum-statistical or hadronic resonance-gas reference. Mean transverse 
momentum per particle can be compared to an inclusive-spectrum reference, and the linear corre- 
lations between these variable pairs can also be compared to an appropriate covariance reference, 
including effects of hadronic correlations. Agreement between data and the GCE might indicate 
that in each event N independent particle momenta are sampled from a fixed inclusive momentum 
distribution (parent), consistent with independent particle production from an invariant source. 
Such a comparison test is closely related to the central limit theorem (CLT) and represents the 
subject of much of this paper. 

4 Equilibration and the Statistical Model 

By the process of equilibration a many-body system in some arbitrary initial state approaches a state 
of maximum symmetry. To study this process with sensitivity one needs a well-defined reference and 
optimized comparison measures. The simplest reference corresponds to an assumption of maximum 
symmetry within a constraint system, which leads to a 'statistical-model' description. Observed 
deviations from a statistical model may then indicate that an equilibration process from a less 
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symmetric initial state has terminated before achievement of maximum symmetry or that the true 
constraint system is incorrectly modeled. 

The specific term 'thermalization' connotes equilibration of energy in contrast to a more general 
situation in which various quantities such as charge, color, flavor, spin and baryon number may 
undergo transport in the approach to maximum symmetry. Equilibration is therefore the preferable 
term. Equilibrium (zero net average transport of conserved quantities) and maximum symmetry 
are essentially synonymous asymptotic conditions. 

4.1 The statistical model 

To establish a robust reference system for detecting possibly small deviations which may convey 
new physics we apply a hypothesis of maximum symmetry (Occam's Razor, maximum entropy): 
we adopt the 'simplest' hypothesis consistent with global constraints and prior information. The 
corresponding many-body system is said to be in statistical equilibrium and its mean- value behavior 
can be described by a statistical model (SM). 

The statistical model has been applied to simple systems including four-momentum and strange- 
ness as additive measures [P, |^, |8|. A dynamical system in some initial state (possibly partonic) 
is partitioned at large scale into elements (clusters, hadronic 'fireballs,' jets, flow elements) having 
arbitrary mean momenta. Small-scale system homogeneity is assumed (all clusters have the same 
intensive properties T and 7^ at hadronization) . The large-scale mean momenta of clusters com- 
prise an unthermalized collective flow field which is assumed to play no role in determining particle 
species abundances. Abundances and ratios are instead assumed to be determined by thermalized 
energy residing at small scale, with no overlap between these two scale regions. 

The partition function for an event composed of individual clusters is a simple product of 
cluster partition functions with no cross terms, implying linear independence of clusters. Each 
cluster partition function in turn treats elements (particles) as members of an uncorrelated gas 
(except for quantum correlations). The treatment is however 'canonical' in that some extensive 
quantities are event-wise conserved, thus introducing a degree of correlation. As we shall see, these 
assumption of homogeneity, linearity and independence are also preconditions for the central limit 
theorem. 

With the exception of the large-scale flow field the system is assumed to be scale-invariant 
and maximally symmetric (minimally correlated). Therefore, the system state is determined only 
by global constraints. Global fitting parameters are T, V and Ag (quark strangeness suppression 
factor) or 7^ (hadronic strangeness saturation factor or fugacity). V is the single extensive quantity 
and should depend on the collision system. T may depend only on the fundamental process of 
hadronization in elementary systems. An asymptotic (maximum) T value may be observed for 
small systems in which the hadronic cascade is quickly arrested. In the limit of large system volume 
and/or particle number this formalism may become consistent with a grand-canonical approach. 

The statistical model provides a satisfactory general description of multiplicity distributions 
(revealing primary-hadron or clan correlations) and hadronic abundances / ratios for small and large 
collision systems over a broad collision-energy interval. An SM treatment of hadron multiplicity 
distributions for light systems |^] results in two-tiered distributions, with primary hadrons following 
a nearly Poisson distribution and secondary hadrons distributed according to a negative binomial 
distribution (NBD) consistent with the number of primary hadrons predicted by the SM. SM 
treatments of hadron abundances or abundance ratios ^ and strangeness suppression Q in 
a variety of collision systems lead to good general agreement with experimental observations in 
terms of a single universal hadronization temperature T ~ 180 MeV and volumes V consistent with 
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specific collision systems from e^-e~ to the heaviest HI collisions. Although the SM provides a good 
general summary of hadronic abundances interesting physics may reside in significant disagreements 
with this model. The SM provides an equilibration reference from which deviations should be the 
subject of further study. 

4.2 Equilibration 

A statistical system appears from a macroscopic viewpoint to move persistently toward a state of 
maximum overall symmetry in an equilibration process. What underlies this apparent teleology is 
understood to be a probabilistic effect: maximally symmetric macroscopic systems are more likely 
to be realized in the random evolution of a normal-system microscopic phase space. Symmetry 
increases with effective phase-space volume (entropy) and dimensionality. Increased correlation 
decreases the effective symmetry, volume and dimensionality of a system and corresponds to a 
less likely macroscopic state. Equilibration then represents a general decrease of system correlation 
toward a maximally symmetric reference as illustrated in Fig. |3|, and can be followed quantitatively 
with differential correlation measures. 

We associate different equilibration types with corresponding transported quantities, such as 
thermal (momentum) and chemical (flavor) equilibration. The relaxation rates of these processes 
depend on local cross sections and flux densities. Different equilibration processes may effectively 
begin and end (freeze out) in different space-time regions of the collision, and with finite transport 
rates the degree of equilibration is limited by the system space-time extension [jl^. An equilibra- 
tion reference should therefore be defined in terms of finite system size/duration to insure proper 
treatment of the scale dependence of system symmetries |jll|]. Concepts such as 'local' and 'global' 
equilibration relate to scale dependence, of vacuum symmetries as well as energy /momentum equi- 
libration. Color deconfinement (spatial color symmetry) may be an established condition within 
cold hadrons. One is more interested in how color and chiral symmetries are modified and spatially 
extended during a HI collision, what is the varying geometry or dimensionality of regions with 
modified symmetries and how does this dimensionality propagate on scale. Microscopic variations 
due to small-scale QCD effects are not as interesting as large-scale deviations from a global thermal 
model that are unique to large collision systems. Comparisons between large- and small-system 
deviations and the scale dependence of correlation within a given system are therefore important. 

4.3 Cascades and scale dependence 

In conventional treatments of an equilibrated many-body system scale invariance is assumed, and 
no provision is made to deal with departures from this assumption. If equilibrium is not achieved 
there has been no quantitative thermodynamic measure system available to deal with deviations. 
The partition systems usually applied represent very elementary versions of Borel measure theory 
and treat scaling issues, if at all, with a limited two- or three-tiered structure. Few-stage partition 
systems in HI collisions represented by jets, strings or primary hadrons combined with a hadronic 
cascade cannot describe more complex departures from equilibrium expected near a phase boundary. 

A classical thermodynamic treatment invokes an extreme scale partition into an 'atomic' scale 
and the system boundary scale. Work can be done on the system boundary 'adiabatically,' meaning 
that energy transferred to the boundary cascades down to the atomic scale within a time small 
compared to the large-scale energy transfer. The details of the cascade are not described. In 
thermodynamic language PdV 'work' represents large-scale energy transfer. TdS 'heat' represents 
small-scale energy transfer. There is no intermediate possibility. What is 'adiabatic' depends on 
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transfer times at large scale relative to transport times (collision rates) at atomic scale in the 
thermodynamic medium. 

In contrast, scale dependence is an essential ingredient of a comprehensive statistical treatment 
of many-body dynamics. A cascade process is a hierarchical structure representing the bidirec- 
tional propagation of correlation and symmetry on scale. Correlation spontaneously propagates 
down scale, and equivalently symmetry (dimensionality) propagates upscale in the dominant trend, 
but with local deviations. In an open system a symmetry increase at larger scale may be imple- 
mented 'opportunistically' by an accompanying symmetry reduction at smaller scale {e.g., turbu- 
lence vortices, biological structures), a process sometimes referred to as 'spontaneous' symmetry 
breaking. 

A falling stone collides with a pool of water both at its surface, depositing energy in surface 
waves, and along a path through its volume, creating a system of vortices. The energy in these 
large- and medium-scale correlations is transported to smaller scale by a cascade process termed 
collectively 'dissipation' or 'damping.' If the dissipation process or cascade is allowed to proceed 
indefinitely the stopped energy of the stone is finally transported to the smallest scale relevant to 
the problem, the molecular scale, and the system becomes maximally symmetric or thermalized. If 
the system is frozen or decoupled before this process is completed residual correlations over some 
scale interval may be observed. The analog in a HI collision is the interplay among momentum 
transfers at nuclear, nucleonic and partonic scales and subsequent evolution to final-state hadron 
momenta at freezeout. 

The evolution of an expanding partonic medium into primary hadrons by a coalescence pro- 
cess 1 12] involves local symmetry reduction (coalescence) accompanying global symmetry increase 
achieved by the expansion. The decay of these primary hadrons to observed secondaries is a further 
cascade process, a symmetry increase in which detectable correlations nevertheless persist in the 
secondary hadron population unless the hadronic cascade can proceed by continued rescattering to 
a state of maximum symmetry - an equilibrated hadron gas. 

The terminology of cascade processes varies with context, as represented by the many terms 
for localized dynamical objects (quasiparticles) in a nuclear collision: fireballs, primary hadrons, 
clans, prehadrons, clusters, participants, (wounded) nucleons, resonances, strings, leading partons 
- resulting in various manifestations of residual correlations in the hadronic final state such as jets, 
minijets, Cronin effect and NBD. It is desirable to develop a simplified common language which 
treats these various manifestations of symmetry variation consistently. 



4.4 Equilibration tests 

Are chemical and thermal equilibrations achieved during a HI collision? Is the final state of an 
N-N, p-A or A-A collision the result of a completed equilibration process? Or are correlation 
remnants from the early stages of the collision process still detectable in the multiparticle final 
state, and do they indicate evidence for QGP formation? Methods are required to distinguish 
between equilibrium and nonequilibrium models. These questions cannot be resolved by inclusive 
measurements. To detect and study possible departures from equilibrium or maximum symmetry we 
require a differential correlation analysis system and a reference system analogous to the statistical 
model but applicable to fluctuation analysis. This paper focuses on measurement and interpretation 
of variance differences as one component of a global EbyE fiuctuation analysis. 

If remnants of a cascade structure are still present after hadronic freezeout we may observe 
the corresponding residual correlations as excess or deficient (co)variance (fiuctuations) of global 
measures. An elementary equilibration test might therefore consist of comparing the variances of 
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event-wise measure-distribution statistics to a reference system. A nonequilibrated system might 
exhibit variances differing from finite-number statistics and elementary hadronic processes. 

Gazdzicki, Leonidov and Roland (GLR) |]l3| considered whether pt production in HI collisions 
is consistent with a fully equilibrated system or retains evidence of initial-state or pre-equilibrium 
scattering of projectile nucleons (ISS, Cronin effect). The initial state scattering (ISS) model pre- 
dicts increased pt fluctuations because pt is generated by a hierarchical process in which transverse 
random walk of projectile-nucleon pt (multiple scattering) is followed by thermal generation of 
hadrons about the resulting primary pt distribution. This relationship between hierarchical struc- 
ture and increased correlation is observed in studies of NBD multiplicity distributions and clans in 
hadron production in small collision systems ISS should lead to a correlated pt population and 
corresponding fluctuation excess. These fluctuations should increase in amplitude with the number 
of possible nucleon interactions (thus with nuclear size). 

Both equilibrium-hydrodynamic and ISS models provide reasonable descriptions of inclusive 
Pt distributions. To further distinguish between the two alternatives GLR proposed to compare 
Pt fluctuations predicted by the ISS model and a hydro model with experimental data. The <I>pj 
measure which they employed is said to provide a quantitative basis for comparing fluctuations in 
HI collisions relative to elementary interactions |jl^. 

An additional source olpt fluctuations according to GLR is the observed N-<pt > correlations in 
N-N collisions |15], some degree of which may survive an A-A collision considered as a superposition 
of elementary N-N events. According to GLR surviving A^-< pt > correlations would lead to 
departures from independent-particle emission detectable with an appropriate measure such as 
^Pt applied in a systematic study of A^-< pt > correlations in N-N, p-A and A-A. They also 
suggest that a substantial change in the nuclear medium in A-A collisions could cause a different 
relationship among elementary N-N collision subevents as indicated by a significant change in 
differential variance. They conclude that should be sensitive both to N-<pt> correlations and 
to ISS. 



5 Fluctuations and Correlations 

Fluctuations and correlations are complementary descriptive terms reflecting the degree of structure 
in a measure distribution. The term fluctuations may carry a connotation of time dependence, but 
the term more generally refers to the variability of a measure distribution, whether temporal or 
spatial. Correlation refers to the localization or size reduction of a measure distribution on an 
embedding space. Increased localization leads to increased variability of the measure, possibly 
interpreted as fluctuations. The terms are sometimes used interchangeably. 

Although fluctuations and correlations can be characterized by a complete system of measure 
distribution moments this paper emphasizes second moments or variances. The primary variance 
source in a typical 'many-body' system is the finite density of measure-carrying objects (degrees 
of freedom or DoF). This concept emerges when a measure system accumulates significant corre- 
lation over restricted scale intervals (correlation onset — > 'particles'), and only a small degree of 
correlation over intervening scale intervals (see Sec. 10. 2| ). As a simplification one then focuses on 



the correlations among DoF at large scale, and ignores their internal structure at smaller scale. 

A reference density of independent DoF corresponds to an expected or basal variance in any 
physically observable quantity. If the effective density of DoF changes from a reference value 
observed fluctuations may depart from expectation. A dramatic change in the number of DoF can 
happen in a phase transition, in which scale-localized correlations are created or destroyed by a 
change in system constraints. 
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Figure 3: Variation of fluc- 
tuations, correlations and sym- 
metry with advancing equilibra- 
tion. Decreasing correlation cor- 
responds to increase of effective 
system size (entropy) within con- 
straints, increased symmetry and 
decreased measure variation (fluc- 
tuations). Thus, determination of 
fluctuations and/or correlations 
of measure distributions informs 
us about the equilibration state of 
a dynamical system. 



Changes in fluctuation amplitude can also happen if correlations among DoF change. Positive 
correlation reduces the effective number of DoF leading to excess variance beyond expectations. 
The opposite is true of anticorrelation. One can say equivalently that the entropy associated with 
a class of DoF is reduced through a coupling process, leading to a degree of freezeout of the 
corresponding DoF from an equipartition process. Changes in fluctuation amplitudes or variances 
can thus be seen as indicators of increased or decreased correlation among DoF or as a change in 
the density of DoF, and can be used to probe the underlying Lagrangian for a many-body system 
or the mechanism of a dynamical process. 

Some correlations in a dynamical system are manifestations of incomplete equilibration, the 
surviving remnants of an incomplete cascade following initial large-scale departures from symmetry. 
Some correlations may be due to quantum statistics or other steady-state multibody interactions. 
If excess correlation has been detected beyond that attributable to known sources the problem then 
lies in identifying the process that led to this correlation and/or what reduced symmetry initiated 
it. From an EbyE viewpoint deviations from expected fluctuation amplitudes (either more or less) 
which may signal unexpected physics or differences from a reference model are of primary interest. 

A central problem for HI physics is the nature of the QCD phase boundary. The location of 
the boundary can be determined by inclusive measurements |l^, but the detailed characteristics 
of the boundary are more directly probed by EbyE analysis [l^. Fundamental to EbyE QCD 
studies is the quantitative measurement of changes in QCD vacuum symmetries as energy and 
baryon-number densities are varied. The most accessible probe of the QCD vacuum may be the 
multiparticle final state of an URHI collision. Near the phase boundary event-by-event dynamical 
fluctuations should provide essential information on vacuum symmetry dynamics. 

In the partonic regime the highly structured nonabelian gluon field plays a substantial role |17|. 
At RHIC energies and above the incidence of parton interactions with high momentum transfer 
should lead to jet formation and a large minijet background in A- A collisions resulting in a 
highly correlated multiparticle final state at high pt due to hard processes. 

At a phase boundary the density of momentum-carrying objects or degrees of freedom (DoF) 
changes rapidly with certain state variables {e.g., temperature and net baryon density). This 
means that in the neighborhood of the boundary a bulk-matter system undergoes a sharp change 
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in pressure and energy density. Fluctuations in extensive quantities exceed those expected from 
finite-number statistics by an amount bounded by (but not determined by) a corresponding linear 
response coefficient {e.g., heat capacity) [0. To determine the detailed properties of the phase 
boundary it is therefore important to search for and measure nonstatistical fluctuations in dynamical 
quantities. 

However, there are important differences between a phase boundary in bulk matter and mani- 
festations of a phase boundary traversed by a small system during a collision, which manifestations 
may be barely detectable in the multiparticle final state [^. Phase-boundary manifestations com- 
pete for observability with details of the collision dynamics, the space-time finiteness of the collision 
system and the varying degrees of equilibration in different space-time regions of the collision. Po- 
tentially observable effects can be attenuated by final-state hadronic interactions which move the 
system further toward equilibrium, with consequent loss of information, and introduce distracting 
hadronic correlations which are not of primary interest. There is a critical need therefore to develop 
a sensitive, complete and interpretable system of correlation measures and comparison measures 
for heavy-ion collisions. 

6 The Central Limit Theorem 

I begin a discussion of variance comparison measures with an elementary treatment of the central 
limit theorem (CLT). In later sections we will develop a more complete view of the CLT and variance 
comparisons applicable to EbyE analysis. The central limit theorem, first stated by deMoivre 
and proved under general assumptions by Markov, plays an important role in many statistical 
applications. It is represented schematically in Fig. ^. The following is an elementary statement 
of the CLT as typically encountered in a statistical context. 

Given a fixed arbitrary parent distribution with mean fi and variance o"^, and given 
an ensemble of events, each composed of A^-fold independent samples from the fixed 
parent: 

• The distribution of event means approaches a gaussian distribution with mean // 
and variance a"^ /N . 

• More generally, the m*'^ moment about the mean of the event-mean distribution 
is reduced from the m^^ moment about the mean of the parent distribution by a 
factor l/N'"'-^. 

The more general {m!'^ moment) aspect of the theorem implies rapid asymptotic approach to a 
gaussian distribution (itself defined by the requirement that all moments beyond the second are 
zero) as — > oo. 

Applied to the problem of event-wise momentum fluctuations in HI collisions as an example the 
CLT can be formulated in the following way. If events formed from N independent p^ samples with 
total momentum Pt (an ensemble of final-state multiparticle distributions resulting from nominally 
identical collision events) are drawn from a static parent distribution (estimated by the inclusive pt 
distribution) the distribution of event means <pf > = P^ /N approximates a gaussian distribution 
with mean approaching the parent mean pi and variance <7<pt> related to the parent variance 0"^^ 
by 



(1) 
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Figure 4: The central limit theorem is illustrated. For an arbitrary parent distribution (left) the 
distribution of A^-fold means approaches with increasing N a gaussian whose width is simply related 
to the parent width. 

Any departures from the CLT are due to 1) (anti)correlated samples and/or 2) variation of 
the parent during the sample history. It is case 2) that we prefer to identify with 'dynamical 
fluctuations,' although this is not a general convention. The term 'dynamical fluctuations' has 
also been taken to represent all fluctuations beyond finite-number statistics |jl^. Net positive 
correlation implies that N^jf < N, that is, the effective sample is less than N. The opposite is 
true for anticorrelated samples. It is important to note that the combined presence of different 
correlation types and dynamical parent variations may nevertheless result in variance behavior 
apparently consistent with the CLT. This ambiguity will be explored below. 



7 The Measure 

To facilitate the quantitative study of fluctuations and correlations in HI collisions Gazdzicki and 
Mrowczyhski (CM) [|l^ developed a differential correlation measure to study equilibration in 
HI collisions. This measure and its experimental applications have stimulated much theoretical 



activity in EbyE analysis recently |20|. CM cite N-<pt> correlations in N-N collisions [^] as a 
possible means to track the equilibration process in heavier collision systems. For the case which 
prompted the development of the measure, if A-A collision events are indeed linear superpositions 
of independent N-N subevents they suggest that a suitably defined fluctuation measure for A-A 
events should have the same value as for N-N events taken singly if linearity prevails. 

According to CM the width comparison measure which they define satisfies two important 
conditions: 1) linearity - superposing independent and equivalent subsystems to form a composite 
system results in the same measure value for subsystems and composite; and 2) zero value for 
independence - if the atomic elements of the subsystems (and thus the composite) are themselves 
statistically independent the measure is zero. This measure has been proposed as a measure of 
EbyE pt fiuctuations. It is also proposed as a general measure to detect event-wise dynamical 



fluctuations in energy and possibly other kinematic or chemical quantities |21|. 

Rather than comparing variances directly is defined in terms of rms quantities with the 
basic structure 



= VNa<p^-ap, (2) 



which is zero if the event ensemble is consistent with the CLT. Following CM one can define an 
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event-wise variation about the mean 



Ze = EC 



Xi X) 
i=l 

= Ne-i<X>e-x) (3) 

and then define the rms comparison measure 

= ^WfN-a^. (4) 

This measure satisfies both of conditions 1) and 2) above. It is zero if the event is composed 
of uncorrelated objects (elements with no resolved structure). The requirement 1), 2) and Eq. 

with <I>a; = are obviously equivalent to the central limit theorem. Given its fundamental 
importance it is not surprising that the CLT is frequently rediscovered. GM apply this general 
comparison measure to single-particle pt distributions as . To specify one replaces x by pt 
in the preceding. 



In a ^>pj analysis of NA49 Pb-Pb data the conclusion was tentatively formed ||2^, 24 1 that 
ISS and N-<pt> fluctuations are substantially reduced, either by collision dynamics and medium 
effects or by final-state rescattering, to a point below statistical significance and certainly below 
values observed in p-A and N-N systems respectively. In a further study ||2^ involving detailed 
simulations of detector response (two-track resolution) and quantum statistical contributions to 
<I>pj a model study was conducted which placed an upper limit on dynamical (parent) fluctuations 
at 1% of the inclusive pt distribution rms variation. 



8 The Structure of 

Application of the CLT in the form of to transverse momentum fluctuations has some unantici- 
pated aspects. As suggested above, there is a number of bipolar and positive-definite contributions 
to If the measure is 'zero' one does not know whether correlations are really absent or instead 
somehow present in equal magnitude but opposite sign from separate sources, or whether bipolar 
effects from a single mechanism distributed over a pair-momentum space make no net contribution 
to this integral measure. Even if is nonzero these ambiguities remain. These issues are not 
confronted in an elementary statement of the CLT as presented above. To better interpret 
and to formulate improved fluctuation measures we examine the structure of in detail. An 
important complicating element is the presence of event-wise multiplicity variations, whereas our 
simple statement of the CLT assumed a fixed sample number A^. One must finally accommodate 
the N-< pt > correlations which originally motivated the development of ^p^. 

Examining the algebraic structure of ^p^ in more detail we will find that this measure is only 
part of a more complex scale-dependent covariance measure system. We will obtain a more compre- 
hensive replacement for the elementary CLT and a generalization of the ^p^ comparison measure. 
This fuller picture will enable us to assess possible physical sources of deviation from the CLT and 
how they are manifest in ^p^ and associated measures. 

Because of potential correlations among statistical quantities and the need for precision in 
obtaining small differential quantities from large numbers with great statistical power I summarize 
the essential elementary quantities. = -M a sum over the event ensemble. J2i = is a 
sum over elements of event e. J2e i I Se = is the mean event multiplicity. Yl,% = is an 
event-wise total measure. X^i^j/Z^i = ^el^e =<x>e is an event-wise mean. X)ei^«/Sei = ^ 
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is an ensemble (inclusive) mean over all events. J2e <^>e / J2e ~ is an ensemble-averaged 

event-wise mean. In general, <x> 7^ x because of possible correlations between N^, and <x>e- In 
fact, 

N(<x> -x) = -r{<x>,N)a<^>aN = cr%c>N (5) 
where the linear correlation coefficient r{x,y) is defined below, and cr%^-^j^ is a covariance. 

8.1 Variances 

This paper emphasizes global fluctuation measures based on variance. The comparison measure 
<I>p^ is based on an rms difference. We now extend the elementary CLT to accomodate event-wise 
multiplicity variation or variation of the sample number N. We consider the interplay among event- 
wise total momentum Pg, particle momentum pi, event-wise mean momentum <p>f. and event 
multiplicity A^e in the spaces Gi{P,N) or Gi{P/N, N). To simplify the notation I omit subscript 
t and indices e and i wherever confusion will not result. 

A differential fluctuation measure based on variances and the CLT can be expressed as 

% = ^(P-pflN-\l{p-pf (6) 



There are two possible ways to expand the variance [P — Pj if the number N of samples or 

particles fluctuates from event to event, depending on how one defines the mean value P. In the 
first case P = Np and I denote the variance by a p. In the second case the mean value P = N^p 
contains the event-wise multiplicity. The variance in the second case is the same Z| appearing 

in the original definition of To simplify I replace (p — p)^ with and Z| with a%. 
The variance of total momentum dp is given by 

al = {P-N -pY 

= Y^j-N{N-l)f + Nal (7) 

The variance used to define the measure is given by 

4 = {p-n ^-pY 

= Y.P'^Pj - - 1)^^ ~ 2p • (NP -N^-p} + Nal 

= al-2p-{NP-N-P}+p'^a% (8) 

This algebra preserves any N-pt correlations. No assumptions are made about the statistical 
interdependence of N and pt- 

8.2 Decomposition of 

We analyze ^p^ as suggested by the structure of (t|. Recalling that ^p^ = ^a'^/N — ap^ we write 

4 = A + B + Nal (9) 
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where 



a%-Nal-fa%, (10) 



and express as 



B = -2p ■ {m <p> - m ■^^} 

= -2p-r{<p>,N'^)-a<p>a,y2 

= -aj, + al+p^a% (11) 



2ap,-N^p, ^ A + B (12) 



This decomposition is nontrivial because of the unique roles played by the resulting A and B terms. 
A is essentially a weighted integral of net two-point correlation , and i? is a measure of N-<pt > 
correlations. Each term is zero under conditions that satisfy the CLT. This decomposition can help 
to disentangle various physical contributions to <I>pj. We now examine A and B in more detail. 

8.3 Two-point correlations 

We provisionally interpret term A as the contribution of two-point correlations to ^p^. In the 
continuum limit this term can be written 



J2piPj-NiN-l)-p^ 

PlP2 ■ {P2{Pl,P2) - Pl{Pl)Pl{P2)]dpidp2 (13) 

The quantity in curly brackets in the integrand is the connected two-point correlation function. It 
is zero in the absence of net two-point correlations. The correlation function used in HBT analysis 
is closely related. The integral is weighted by the prefactor piP2^ which leads to an undesirable 
emphasis on the large-p region. Contributions to A include, but are not limited to, quantum 
statistics. Coulomb effects, resonance decays, instrumental effects and 'dynamical' fluctuations. 

It is more useful to study the two-point correlations in Eq. (|l^ directly rather than in terms 
of an integral weighted by the factor pip2- One can switch to a mixed-pairs reference 

Pl{Pl)Pl{P2) Pmixed{pi,P2), 

and compare the two-point densities on {pi,P2) for sibling pairs and mixed pairs. A then becomes 

A ^ J piP2 ■ {psib - Pmixed}dpidp2 (14) 

We can plot the ratio of density distributions Psib/ Pmixed for sibling and mixed pairs, and also 
project the raio onto momentum difference q and mean momentum k by analogy with standard 
HBT analysis. 

To minimize instrumental effects and to reveal possibly different correlation mechanism for 
different pair types we plot the pair density ratio Psib/ Pmixed in Fig- H |26| separately for two 



combinations of like-sign pairs, one of unlike-sign pairs and one for all pairs taken together. One 
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Figure 5: Distribution of sibling 
to mixed pair-density ratio on pt<^ 
Pt space for 158 GeV/nucleon Pb- 
Pb collisions. Density-ratio vari- 
ation (denoted by box size) is 
within the interval 1 ±0.003. The 
histograms are symmetric about 
the diagonal by construction. The 
pair momentum difference is q 
and the momentum sum is k. The 
variables pi are transformed from 
rrit — mo in such a way as to 
make the statistical error approx- 
imately uniform over the two- 
point space. 



observes a prominent ridge along the main diagonal {q = 0) for the like-signed pairs corresponding to 
quantum correlations and a peak in the unlike-signed pairs (k ~ 0) due to the Coulomb interaction. 
There is also a very broad aniicorrelation at large q in the like-signed pair distributions and a broad 
correlation at large q in the unlike-signed pair distribution. These last results are unexpected and 
illustrate the difficulty of interpreting in a simple manner. A separate treatment of positive, 
negative and unlike-signed pairs is generally necessary in order to extract complete information. 
Different physics may affect each combination. 

Evidence for dynamical fluctuations in the parent distribution should be common to all pair 
types (like and unlike pairs), providing an additional signature of true dynamical fluctuations. 
Because is an integral correlation measure there could be major differences between the psih and 
Pmixed density distributions which nevertheless remain undetected using Scale-local measures 
and direct analysis of individual pair spaces should provide a much more sensitive approach to the 
search for dynamical fluctuations. 

8.4 N-<pt> correlations 

Term B is proportional to a covariance or linear correlation coefficient. 

B = -2p ■{N'^ <p>-W- <p>} 
= -2pr{<p>,N^)Na^>aN 

= -2p^^>Ar2 (15) 
where the linear correlation coefflcient (Ice) for two variables x and y is defined by 

{x-x){y-y) 
r(x.y) = 
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and the numerator is the covariance a^y. 



(16) 




Figure 6: Illustration of A'^-< pt > correlations in light systems (left) and HI collisions (right). 
N-<pt> correlations observed in light systems (N-N) may be used as a probe in HI collisions to 
study equilibration. The question then becomes to what extent N-<pt > correlations survive the 
equilibration process; how are conservation laws expressed differently in HI collisions than in N-N 
collisions? 

Note that the Icc also has properties 1) and 2) used in the derivation of $p^. It is 1) invariant 
under linear composition of independent and equivalent subevents and 2) zero in the absence of 
correlation of event elements. To demonstrate the linearity of this measure consider Pb-Pb collisions 
as linear superpositions of N-N collisions. If fcg N-N subevents are independently superposed in 
Pb-Pb event e, < Pt >e is the event-wise mean pt, <pt >j is the mean pt for the j*'* N-N subevent 
and rij is its multiplicity then: 



r{<P>,N)pb_Pb = 



i<Pt>-Pt)-{N-N) 

cr<Pi> ■ <^N 
E-,j=i <Pt>i -rij-TfN 

k ■ <pt > ■n + k{k — 1) ■ptn — k -pin 

= r{<p>,n)N_N 

Thus, the Icc also offers a significant opportunity to probe 1) momentum conservation effects, 2) 
equilibration processes and 3) independence of particle emission. 

The Icc is linear under composition of partition elements (scale invariant) because it is a ratio 
of total covariance matrix elements - a'^y/axCTy A scale-dependent factor M cancels in this ratio. 
If the system is scale invariant (satisfying CLT conditions) then this Icc ratio is invariant also. 
However, if the CLT is not satisfied this composite measure is not optimal for understanding cor- 
relation sources, because it mixes the scale dependence of variances and covariances. It is better to 
look at cr^^j^, (7^ and cr^ or their total-covariance counterparts separately. The intuitively motivated 
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linearity condition for will become a generalized scale-invariance test of total-covariance ma- 
trix elements (cf Sec. 10.2 ). The question of the evolution of N-N N-<pt> correlations in Pb-Pb 
collisions which was raised in [|l^] and [|l^] is characterized in Fig. 



9 A Vector Model of Variance 

We find that the variance comparison measure defined in the spirit of the CLT apparently has 
two algebraic components, A being the integral of net two-point correlation and B containing the 
Icc for N-<pt> correlations. In summary 

al-Nal = A + B 

« 2Nap<^p, (17) 

We require a clearer picture (a geometrical model) of how the quantities Up, (t|, and aj^ are related 
to one another and to the CLT. This problem involving measures P, N and P/N is a particular 
instance of the generic case (mi,m2, < m >= mi/m2). 

I now develop a vector representation that compactly displays the algebraic relationships among 
the variances. In a later section I adopt a broader description based on the covariance matrix. To 
simplify notation I use a system of commensurate variables {P,pN, N <p>), based on the original 
{P,N, <p>) but all having mean value P and respective rms widths {paN,crp, Na<p^ = (t$). A 
sketch of typical 2-D density distributions on these variables with corresponding variances is shown 
in Fig. 1^ For discussion I label these single-point measure spaces Gp and The distributions 
on Gp and are simply related algebraically. For events consisting of independent samples the 
covariance in G$ is small and that in Gp is therefore maximal. As represented by the dashed ellipses 
in Fig. 1^ the opposite could also be true in principle, or any intermediate case could occur. This 
descriptive problem underlies the $pj-variable analysis. I now develop a vector model of variance 
which can be used to relate the variances and covariances in these spaces. 



9.1 A vector covariance representation 

Because of their simplicity the density distributions typically encountered in spaces Gp and G,^ 
for an approximately homogeneous system can be well described in terms of their first and second 
moments (gaussian fiuctuation model), with the addition of a covariance element describing the 
correlation between fluctuations on two degrees of freedom. The covariance matrix in either space 
can be written formally as a direct product of column/row matrices whose elements are themselves 
vector representations of variance defined in a pseudospace. For example, we can define a generic 
covariance matrix 
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where the dot product dx • dy is defined by 

dx ■ dy = Ox (TyCOSO (19) 
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Figure 7: The measure space on the left is Gp = Gi{p,n) and on the right is G$ = Gi{p/n,n). 
These spaces serve as archetypes for any related set of of extensive and intensive (ratio) measures 
or thermodynamic state variables. The solid ellipses illustrate the case of independent sampling at 
small scale followed by composition to events, whereas the dashed ellipses illustrate independent 
sampling at large scale followed by partition to N objects. The major axes of the ellipses include 
variations of the effective event volume, while the minor axes are determined primarily by statistical 
variance. 



and cos6 = r{x,y) is the linear correlation coefficient between variables x and y. 

This system can be used to relate the variances in spaces Gp and To do this I first form a 
linearized relationship between the two spaces. I then have 

N{P/N -p) PS {P -P) -p{N -N) (20) 

with corresponding variance vectors related by 

(T$ « ap—puN (21) 

(or more generally ?Ti2 i?<m> = (Tmi — fhi/fh2 ■ ffmz)- This vector relation is illustrated on the 
left side of Fig. |8|. Examples of global intensive quantities for a hadronic system are Nk/Nj^ 
(flavor fluctuations), Nj^/N^ (isovector fluctuations), iVo/(iV+ + A^_) (isotensor fluctuations) and 
N^_§/Nt^ (baryon-number fluctuations). 

With some manipulation of Eq. ( pi] ) we obtain two useful relationships 

(Tp = cr| + 2ri paN 0"$ + p'^ajf (22) 

(t| = CTp - 2r2 paN crp + p'^a% (23) 

where ri = r{<p>, N'^) and r2 = r{P,N) are the linear correlation coefficients. I define angles 9i 
and 02 by ri = cosOi and r2 = 00362. The geometrical interpretation of these angles is also shown 
in Fig. |8|. By projecting out cartesian components I obtain useful sine and cosine equations 

paN = ap 00362 — (T^ cosOi (24) 



(T$ 3in6i = ap sin92 = ViV dp + 5a 



(25) 
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Figure 8: A vector representation 
relating variances in Gp and G$ is 
shown above. aR represents a refer- 
ence, and fjp, £7$ equivalently rep- 
resent the data. The dependence of 
the two variance quantities on a hn- 
ear correlation coefficient is shown at 
right. 
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the second of which includes a possible contribution ba from an unspecified variance source. Eqs. 
( p^ and (24) are analogs to Eqs. (17) and (23) respectively of [^]. In the present case the results 
assume no special conditions and arbitrary deviation from the CLT whereas the results in [16| are 
derived in the specific context of an equilibrated pion gas. 



9.2 Vector deviation from a CLT reference 

The vector description of variance in spaces Gp and has so far been generic. That is, —cosQ\ 
and cos^2 could have taken any consistent pair of values on the interval [0, 1]. The two spaces are 
complementary and interchangeable. A reference system for the CLT applicable to HI collisions can 
now be defined by fi/j = [per at, \/iVcTp]. This choice implies zero covariance between variables <pt> 
and [cosOi = 0) and momentum variance given by a CLT expectation The number variance 
o"^ exceeds the Poisson expectation N primarily due to correlations from resonance decays. This 
reference is consistent with the central limit theorem but represents a more complete description 
(including covariance). The physical interpretation of this reference is consistent with the CLT - 
individual particle momenta chosen independently from a fixed parent distribution represented by 
an inclusive pt distribution. 

Given this reference choice I can define a vector measure of deviation from the CLT reference, 
vector Act running from the reference (ct_r) to a point representing the data (cxp) 

Act = Bp — ap 

= [a^cos9i,a^sin6i — \fNap] 

^ VN[ap 591,5(7] (26) 
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where we finally retain only terms linear in small quantities {59i,6a). We see that two quantities, 
69i ~ cosOi = Icc oc B and 6(j oc which both have the properties 1) - linearity and 2) - zero for 
independent objects, provide a measure of deviation from a dynamical variance reference relevant 
to the CLT. This extends the initial concept of the measure finally to handle also A^-< pt > 
correlations. 

Dependence of the variances on ri = cos9i is shown on the right side of Fig. ^. Note that 
depends quadratically on r{<p>, N^) = cosOi, which is just the measure sensitive to N-<p> 
correlations. This dependence was first suggested by numerical simulations and implies that 
<&P( itself is ironically quite insensitive to N-<p> correlations. Combining Eqs. ([TtI), (^) and ( p5|) 
I have 

= a,;>/VN-ap 

~ ^cos^0^ + Sa. (27) 



which is consistent with Eq. (2£). We conclude that for nearly symmetric systems {ap ~ aji) 
is primarily sensitive to variance sources (dcr), whereas covariance sources such as N-<p> 
correlations must be measured explicitly by a covariance or linear correlation coefficient. In Sec. 



11 we will consider the relationship between covariance matrices in Gp and G$ and similar spaces 



more generally. 

A concern is raised in |24] that cr^^> = ^Pt/N ^1/-^ ' that is <I>pj is not directly related to 
(^%p£^- But it is easy to show that the latter quantities in the present case are arbitrarily close, 
to the extent that the relationship between spaces Gp and G$ can be linearized. This depends in 
turn on the degree to which the CLT conditions are met and on the size of the total multiplicity 
N (equivalently the scale interval from particle to event). 



9.3 Significance of A, B and $ 



Pt 



We now reexamine A, B and in light of this more complete vector variance representation. 
The relationships among A, B and were first suggested by numerical simulations [27| which 



demonstrated a strong correlation between A and B when N-<pt > correlations were artificially 
introduced into a mixed-event reference. A and —B varied together over a 40MeV/c range while 
<I>pj fluctuated about zero. This suggested either that A and B were redundant or that one carried a 
subset of the information contained in the other. The relationship A + B 2Nap^ simply defines a 
plane in a 3-D space of {A, B, <I>) without giving emphasis to any one of the three related quantities. 
We would like to develop a more substantial understanding of this relationship. 
Given the definitions of A, B and <I> I can express A as 

A (t| — Nap + 2paN Vn ap cos9\ 

2VNap- \/N^ + 2paN ■ VNapCosOi 
w 2Aa ■ ap (28) 

A thus has the form of a vector dot product. Note that if the deviation Aa is orthogonal to the 
reference ap there is no net two-point correlation, even though the CLT is not satisfied. We now 
have a relationship between departures from the CLT on the one hand [cosOi and ^p^), which are 
parts of a global measure system, and net two-point correlations on the other {A), which refer to the 
microscopic particle complement. This vector result suggests that A is in a measure class distinct 
from cosBi and $pj. We will see below that this is indeed the case. In pursuing the structure of 
$P( we have discovered a broader organization. 
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We also find that the presence of two-point correlations may result in deviation from a refer- 
ence variance system, but not necessarily. On the other hand deviation from a reference system 
may result in net two-point correlation, but again there is no necessary relationship. Variance 
changes and N-<pt > correlations (covariances) are global effects with possible thermodynamical 
interpretations. Both may correspond to net two-point correlations at the microscopic level or vice 
versa, but there is no necessary connection. This picture is still not fundamental. We require a 
more complete treatment of variance and correlation, including explicit scale dependence, to fully 
interpret any results. 

This treatment of variance provides a geometrical model relating archetypal spaces Gp and 
It illustrates how small incremental contributions to variance and covariance relative to a reference 
are related to the measure. It illustrates dramatically the reason that is insensitive to A^— < 
Pt > correlations, which are instead measured by a covariance element. We have established that 
the decomposition of into components A + B is nontrivial, that in fact itself approximates a 
component of a more comprehensive measure system, and that its value is not related in a necessary 
way to the presence or absence of system correlations, either differential or integral. To proceed 
further we must extend our description to a system of scale-dependent covariance matrices. 

10 Scale-Dependent Variance 

To address the ambiguities encountered in the ^p^ analysis we now attempt to extract more com- 
plete information from collision events. We do this by expanding the description in two ways: 1) 
we move to an explicit scale-dependent description, and 2) we adopt a description based on the 
complete covariance matrix for all additive measures relevant to the dynamical system. The vec- 
tor representation of variance developed previously is contained in a larger description involving 
scale-dependent covariance matrices and correlation integrals. 

There are four spaces involved in a variance study. -Pi(x) is the single-point primary space 
with measure-density system fh defined on it. Q2{x®x) is the corresponding two-point space with 
measure product density m®fh defined on it. These spaces may not be directly accessible to obser- 
vation. Measure integrals for specific partition systems are themselves distributed on spaces Gi{rh) 
and G2{ifn m) which may represent the results of direct measurement. The covariance matrices 
described below represent a gaussian approximation to distributions on Gi{rh). Covariance- matrix 
differences are in turn related to correlation integrals AC2(?7t-i, m2; (5x2) over two-point space 
Q2{x ® x) and covariance integrals A{5xi,5x2'-, mi, m2) over space G2{m m). 

10.1 Variance and correlation integrals 

A variance is implicitly defined at a particular scale (the bin size of a binning system). Each entry 
in Fig. 1^ corresponds to an event, which can also be considered as a partition element or 'bin' in 
some larger space (Sec. ^) resulting from the sampling of a general thermodynamic system at some 
arbitrary scale. We now move to a more explicit and general scale-dependent binning approach, 
with events and particles considered as particular cases. 

I assume a space in which particles, events and event ensembles are arrayed, ft is the primary 
or P space. This space is binned at some arbitrary scale 5x. At some smaller scale a bins may be 
identified with particles. At a larger scale L bins may coincide with events. The space is arbitrarily 
large and may contain many events (an event ensemble). Associated with any P space is a set of 
Qq spaces consisting of (7-fold cartesian products of P space. In the discussion of variances that 
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follows the space Q2 will be of primary importance. This space represents all the point pairs of P 
space. 

I consider relationships among measures on a bounded scale interval [(5xi, 6x2] on which I define 
a scale-dependent covariance matrix K(m,(5x). (m is the set of measures defined on fi). This 
scale-local generalization with particles and events as limiting cases provides a more complete and 
self-consistent description of variance. I begin with the scale dependence of variance for a single 
measure. 

Consider an additive measure (say momentum) distributed on and binned at two scales 
{5xi,5x2), with typical bin contents p{5xi) and p{5x2) and with p{5x2) = J2iLiPii^xi) for the 
smaller bins contained in a larger bin. I can then write for each large-scale bin 

N N 

If I average over all bins I have 



p'^{6x2) — N p'^{dxi) = ^'^pi{6xi)pj{6xi) (30) 

where N{6xi,6x2) = (6x2/6x1)'^ = M{5xi) /M{5x2) is the number of smaller bins contained in a 
larger bin {d is the P-space dimension, which I hereafter omit to simplify notation). The RHS is 
also expressible as the integral over pair-momentum space of a two-point density distribution. I 
discuss this integral in more detail below. The LHS can be expressed as the difference between two 
rank-2 correlation integrals. 

The rank-g normalized correlation integral at scale 6x for an additive measure p can be approx- 
imated by 

M{5x) ( M{5x) 1 ^ 

C,{5x) ^ P'(^^)/ E P^(^^) (31) 

i=l { i=l J 

where M(dx) is the number of occupied bins at scale 5x in the space Q (containing nonzero measure). 
This bin-based definition approximates an integral on scale in space Qq from zero scale up to the 
bin size 5x. For an uncorrelated or uniform distribution Cq{5x) Cq,ref = M^~''(5x), which is 
one possible reference choice. 

I can express mean-square quantities in terms of the correlation integral by 



p^{5x) « M{5x)-f{5x)C2{5x) (32) 



from which I obtain 



p^{5x2) - Np^{5xi) ^ M{6x2)p\6x2)C2iSx2) - [M\Sxi)/M{6x2)]p\5xi)C2{5xi) 

« P\6X2) M{5X2){C2{6X2) - C2{5xi)} 

« p\6x2)M {6x2) 02(6x1,6x2) (33) 

The difference between correlation integrals at two scales is just a correlation integral over a scale 
interval, 02(6x1 , 6x2)- This justifies the observation made earlier that at least some global variables, 
those based on moments, are scale integrals. If we subtract from both sides of Eq. ( |3^ ) the 
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corresponding expressions for a reference distribution such as a uniform distribution, with 

as shown in Fig. 



P 



and Y.ijLjPiPj = N{N 



then we have the variance difference 
p'^{5x2)M{5x2) AC2{p,p,6xi,6x2) 



(34) 



where we adopt the notation (mi, 1712, dx) to anticipate generahzation to covariance matrices for 
arbitrary measure systems. AC2{p,p', 6x1,6x2) is then the difference between correlation integrals 
over a bounded scale interval for object and reference distributions. The implicit reference for a 
standard variance is a uniform distribution on the system of occupied bins or support of the object 
distribution. The variance quantities a'^{p,p,6x) could also be defined in terms of an arbitrary 
nonuniform reference. 




Figure 9: Two-point space Q2{x 
x) corresponding to a 1-D primary 
space. Statistical quantities are 
given for a uniform reference distri- 
bution. Q2 is binned at three arbi- 
trary scales. A comparison is made 
in the text between the variance con- 
tained in bins at a smaller scale 
(dark) and at a larger scale (light). 



Fig. ^ illustrates several of these concepts. Bins at smaller (dark) and larger (light) scale are 
highlighted along the main diagonal of Q2{x ® x) space. The variation of total variance (next 
section) with scale implicit in Eq. (p4|) can be illustrated with the following argument. a'^{p,p, 6x1) 
represents the variance per bin for the smaller bins, while a'^{p,p, 6x2) represents the same quantity 
for the larger bins. Since there are smaller bins in a larger bin a'^(p,p,6x2) — Na'^{p,p,6xi) 
represents the excess variance in a larger bin not accounted for by that contained in the smaller 
bins. This excess must be contributed by the N{N— 1) off-diagonal smaller bins contained in a larger 
bin. Within a factor, the measure of additional variance in the off-diagonal bins is the correlation 
integral difference /S.C2{p,p]6xi,6x2) over the scale interval bounded by the two bin sizes. The 
two diagonal arrays of bins approximate the integration regions (diagonal strips) appropriate for 
exact correlation integrals C2{6xi) and 02(6x2) in Q2. This connection of variance with correlation 
integrals is in part motivated by subsequent integration of covariance treatments with scale-local 
topological measures P, p^ ]. 

I now return to the cross term 



P\P2 ■ P2{P1,P2) dpidp2 



(35) 



This relation connects the two-point measure space G2{p,p) (right side) with the two-point space 
Q2{x,x) represented in Fig. ^ (left side). The integral over the entire two-point momentum space 
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of the density P2{Sxi,Sx2;pi,P2) is equal to a correlation integral over a specific scale interval in Q2 
of the density P2ip,P] xi, X2). Distribution p2iSxi, 6x2] Pi, P2) is a parametric function of that scale 
interval as indicated by the argument. Note the convention p2{params; variables) where params 
determine the nature of the density and variables are the position indices for the corresponding 
space. 

It is useful to adopt the convention that symbol a represents a 'particle' scale and L rep- 
resents an 'event' scale, with il. representing the entire containing space. In the particle-event 
limit N{6xi, 6x2) — > Npart, and M(6x2) — > M. (total event number) as (6x1,6x2) {a,L). The 
distribution in Gp = Gi{p,n) is then p2{L,^};p,n), and the two-point particle momentum distri- 
bution in ^2(^1,^2) is p2{a, L;pi,p2)- We can generally write p2{6xi,6x2]Tni,mj) for the density 
in G2{mi,m2), consistent with the notation above. To complete the treatment we write 

a^{p,p,6x2) - Na'^{p,p,6xi) = 




P1P2 ■ {P2,obj{Sxi,6x2;) - P2,ref{5xi,6x2])]dpidp2 

= A{6xi,5x2]p,p) (36) 

which defines A{6xi,6x2'-,p-,p) as an integral over the two-point measure space G2{pi,P2)i the same 
A which appears in the decomposition of $pj. In case p2^refi^xi, 6x2;) corresponds to a uniform 
reference then A{6xi,6x2]p,p) is the covariance of the two-point density p2,obj{6xi, 6x2;pi,P2)- 

This result for variance a'^{p,p, 6x) applies also to other elements of any covariance matrix 
for a measure system m. Thus, the covariance difference a'^{p, 11,6x2) — Na'^{p,n,6xi) equals an 
integral over the space G2{p, n) of a corresponding two-point density difference P2,obj{^xi,6x2;p, n) — 
P2,ref{6xi, 6x2;p,n) represented by A{6xi,6x2;p,n). 

10.2 Total variance 

The linearity property invoked in the derivation of ^pt is actually an implicit invocation of scale 
invariance. The approximately invariant measure is the total variance introduced in the previous 
section. Recalling that N{6xi,6x2) = M{6xi)/M{6x2) and multiplying through Eq. ( |3^ by 
M{6x2) (omiting subscripts) I have 

i:'^{p,p,6x2) - 'S'^{p,p,6xi) PS AC2{p,p;6xi,6x2) (37) 

where T,'^{p,p,6x) = M{6x) a'^{p,p,6x) is the total variance of a measure p on Q relative to a 
uniform reference at scale 6x {a'^{p,p, 6x) is the per-hin variance), and P = M{6x)p{6x) is a scale- 
invariant total measure on Vl. We establish here the concept that variation of total (co)variance 
on scale is the underlying issue for the central limit theorem and any correlation analysis based on 
global variables. Variation of the scale-dependence of total variance may be due to changes in the 
underlying system Lagrangian for certain constraint configurations (the neighborhood of a phase 
boundary), in which a correlation onset may appear or disappear, producing large fluctuations in 
the primary DoF which determine the basic statistical fluctuations for all observables. 

We can connect increased total variance with its source as represented by the correlation integral 
by developing a differential equation. We simply take Eq. (37) to the limit of small scale interval 



dT?{p,p,6x)/d6x « P'^ AdC2ip,p;6x)/d6x 



(38) 
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Figure 10: Scale-dependence of to- 
tal variance for a system with scale- 
localized correlation onsets. Correla- 
tion onset extends from the charac- 
teristic size of an object to the mean 
spacing between objects where local- 
ized objects predominate. The dis- 
tribution of total variance is not nec- 
essarily monotone increasing with 
scale due to the possibility of an- 
ticorrelation in distributions. Such 
curves may vary with global con- 
straints such as (T, fi) near a phase 
boundary, leading to departure from 
a CLT expectation. 



where dC2{p,P', 5x)/d5x is an autocorrelation density. This relation then implies that any variation 
in total variance over a scale interval (the object of a CLT test) corresponds to a nonzero difference 
between object and reference autocorrelation densities over that interval. Since the autocorrela- 
tion density is a quantitative symmetry measure we see that scale-invariance of total variance is 
intimately related to the scale-dependent relative symmetry of object and reference distributions. 

The CLT provides a useful reference for measure systems in which localized regions of correlation 
onset are separated by significant scale intervals with negligible accumulation of total variance 
(approximate scale invariance) as illustrated in Fig. Distinguishable objects that result from a 
localized correlation onset can be designated 'particles' or dynamical DoF, and correlation analysis 
based on integral quantities (global variables) is then simplified to charcterizing any (relatively 
small) total variance accumulated over some scale interval between these correlation onsets. 

Approximate invariance of total variance over a bounded scale interval is espressed by T,'^{5x2) ~ 
E^(fei) for some [6x1,6x2]- We then have 

N{6xi,6x2)(T%p^^ = crpjN{6xi,6x2) 
= T?{5x2)/M{5xi) 
« T.^{5xi)/M{6xi) 

which is a statement of the CLT as used in the definition of the variable (Eq. (|^)). 

In summary, the accumulation of total variance on scale is a basic issue for any global- variables 
analysis based on a CLT test. This is simplest to study in cases where one has scale- localized 
correlation onsets separated by scale intervals of nearly constant total variance, over which a CLT 
reference is a good approximation. We have shown - Eq. (|3^ ) - that the difference ap — N cFp is 
actually a scale integral represented by A[6xi,6x2]p,p)-, and that it is desirable to examine the 
integrand of A directly in a differential analysis whenever possible. An extention of this program to 
scale-local topological measures examines the detailed scale-depence of total- variance accumulation 



(39) 
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or its equivalent without prior assumptions With a scale-dependent treatment of variance in 
hand we now proceed to a general system of scale-dependent covariance matrices. 

11 Scale-Dependent Covariance Matrices 

We now consider an arbitrary system of dynamical measures and the corresponding scale-dependent 
covariance matrix. A covariance matrix provides a gaussian approximation to a scale-dependent 
density distribution in measure space Gi{m). As we have shown, variation of (co) variances on 
scale can be equated to corresponding integrals over two-point measure spaces. We distinguish 
between correlation integrals of densities /02("ii; "12; xi, X2) on spaces Q2{xi,X2) represented by 
C2{mi,m2] 5xi, 8x2) and integrals of densities p2{5xi, 6x2; nii, 1712) on two-point measure spaces 
G2{mi,m2) represented by A{6xi,5x2;mi,m2)- 

11.1 Single-point and two-point measure spaces 

If a space Pi{x) supports a set of additive measures m = (mi, ?Ti2, ...) (e.g., momentum p and 
multiplicity n) having bin values mi{5x) for the ith bin of size 5x represented as a density in a 
single-point measure space 6*1(771) with projections onto pairs of measure values such as Gi{p,n), 
and I form a two-point space Q2ix0x) (Fig. 1) then I can define a corresponding general two- 
point measure space G2(mi 070,2) with projections onto pair spaces such as G2{pi,P2), 62(^1,712), 
G2ini,n2) and so forth. Each bin in Q2 or pair of bins in P space contributes one 'point' (integral of 
the measure densities over the bin) to a density distribution p2iSxi, 5x2; mi, 7712) for a specific scale 
interval in each of the two-point measure spaces (these spaces are binned also, at the resolution of 
a detection device) . These density distributions are parameterized by the scale intervals in P space 
used to define them. 

Specifying a scale interval {5xi,5x2) may seem mysterious. The lower limit is obviously the bin 
size in P space used to define the densities. But to form a correlation space one must specify over 
what region of P space bin contents will be correlated (Fig. ^). The size of this region (a type of 
correlation length) provides the upper limit to the scale interval. In an EbyE analysis 'particles' 
are correlated within 'events' and not between events over an entire ensemble. 

Measures defined at the event scale may be used to sort events into categories. One aspect 
of EbyE analysis is to establish such categories. Subsequent EbyE analysis may indicate that a 
measure space is incomplete by revealing previously unknown correlation structures within events. 
In defining references for correlation analysis it is desirable to form reference pairs from events that 
are 'nearest neighbors' in this expanded measure space, for instance, events with the closest total 
multiplicities or transverse energies. 

The quantity A(5xi,6x2',p,p) defined in Eq. (^) is an example of an integral over a two- 
point measure space. It is just the term A identified in the decomposition of the <I>pj measure. 
It is the integral over measure space G2(pi,P2) of the connected two-point correlation function 
P2,obj{(i, L;pi,p2) — P2,ref{ci', L;pi,p2) (Fig. |5|). Similar integrals for various additive measures taken 
pair-wise form the elements of the integal matrix A{5xi,5x2]m ® m). The relationship between 
distribution moments and g-point densities is also discussed in [29|. 



11.2 Covariance-matrix elements 



Given a general measure system fh defined on and arguing by analogy with the momentum- 
variance results above, I can write for any pair of component measures mi and m2 the quadratic 
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mean 

M {Sx)mi{6x)m2{Sx) = M'^mifh2C2{'rni,m2;Sx) (40) 

where 02(1711,1712; Sx) is the correlation integral up to scale 6x of the two-point density distribution 
P2ii7ii,m2] xi, X2) on Q2- Note that on the left side (mi, 771,2) are measure values from the same 
bin in Pi whereas on the right they come from different bins (off-diagonal element of Q2)- The 
corresponding total (co)variance is given by 

M{6x)a^{mi,7n2,Sx) = M{6x) {mi{6x) m2{dx) — mi{5x) ■ m2{dx)} 

= T,'^{7ni,7n2,Sx) 

= M'^7hifh2AC2{mi,m2; 5x) (41) 

A comparison of quadratic means at different scales is given by 

M{6x2) 7ni{6x2)7n2{Sx2) — M{6xi) mi{6xi)m2{5xi) = M^mi7fi2C2(mi, 7712; fei, 5x2) (42) 

and a (co)variance comparison (generalized form of CLT) is given by 

o"^(7n,i, 7n,2, 6x2) — N {6x1, 6x2)cr'^ {7171,1712, dxi) = rfii ■ 7712 M{5x2) AC2(7n,i, 77^2; 5xi, 8x2) 

= J 7ni7n2{p2,obji8xi,6x2;mi,m2) - p2,refiSxi,6x2;mi,7n2)}dmid7n2 

= A{6xi,6x2;7ni,m2), (43) 

or in terms of total (co)variance by 

(7711, 7772, 5X2) — (7771, 7772, fei) = M^?fli777,2 AC2 (777i , 7n,2 ; , (^^2) 

= ^((5x1,(5x2; 7771, 7712). (44) 

The quantities (7771, 7772, 5x) are elements of a total covariance matrix IC{m,5x), The per-bin 
variances and covariances tr^ (7771, 7712, (5x) are elements of the usual covariance matrix K(7n,, (5x). 
The quantities 777,1 • 7772 M((5x2) AC2 (7771, 7772; (5xi, (5x2) are equal to corresponding elements of the 
integral matrix A{5xi, (5x2 ; ttt, ® 777) , and the quantities M^777 i777,2AC2(777i, 7772; 5xi, 6x2) are equal to 
elements of an integral matrix A{5xi, 6x2; m ® 777), both corresponding to a specific scale interval. 
We can summarize with the matrix equations 

K(777, (5x2) — A^((5xi, (5x2) K(777, (5xi) = A((5xi, 5x2; 777 777) (45) 

and 

/C(?77, (5X2) — /C(777, (5xi) = ^((5xi , (5x2 ; 777 777,) (46) 

which relate covariance matrix elements and integrals over two-point measure spaces for a bounded 
scale interval. 

11.3 Covariance matrices in Gp and 

We now apply these new concepts to the original problem of understanding the relationship among 
A, B and ^p^. The spaces Gp = Gi{p,pn) and G<i> = Gi{np/7i,pn) « Gi{p — pn,pn) shown in 
Fig. ^ provide an archetypal system for the study of extensive/intensive variable combinations 
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and associated covariance structure. These spaces are approximately related by a linear shear 
transformation S. The corresponding covariance matrices Kp and K$ for an arbitrary distribution 
common to the two spaces are then defined by the relations 

= 




(47) 

These matrix equations compactly represent the relationships among covariance matrix elements 
derived in Sec. using a vector representation. 

I now define a reference covariance matrix at small (particle) scale in G$ which is analogous 
to (JR in the vector variance formulation and which may serve as a null hypothesis based on a 
physical model. I obtain the corresponding form in Gp by applying the same shear transformation 
as before. The reference matrix in G$ has zero covariance and variances ap and o"^ determined by 
some unspecified elementary phenomena. For 6xi ~ a, ap may correspond to a thermal distribution 
and cr^ ~ 1 modulo quantum statistics. The momentum p represents a characteristic momentum 
at small scale which may coincide with p. 

K,(*..) = ( f ) (48) 

This choice is particularly simple and could represent an independent-particle-emission hypothesis. 
We can now obtain expressions for covariance difference matrices in both spaces. 

Using the inverse transform K.p(dx) = S^^ • K$((5x) • S""*^ I can write 

,2 Ar_2 -2Ar_2 -2 ^2 at ^2 \ 



( 4 - + 2j5a^>^2 pa%p^^^+p^{al- Nal) 



(50) 



I can now relate covariance matrix elements in the space Gp = Gi{p, n) to corresponding two-point 
measure integrals in the space G2{{p,n) (8) {p,n)) by the matrix equation 



K{p,n,6x2) - N{6xi,5x2)K{p,n,5xi) = A{5xi,5x2; {p,n) (Si {p,n)) 

( -^pp P-^pn \ 



P-^pn P^Aj-, 



(51) 



The LHS is the scale difference between covariance matrices for two-point densities defined on 
subspaces of Gi{m) (measure values taken from the same bin). The RHS represents covariance 
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integrals for two-point densities defined on subspaces of G2{rh®fn) (measure values taken from dif- 
ferent bins), the most integral correlation measures that can be defined on G2. Thus, (co)variances 
of measures defined on single bins in a primary space are interpreted by two-point correlations 
among measures on the primary space as determined by covariances on pairs of bins. 

We now return to the original problem of A, B and <I>pj. Combining Eqs. (^) and ( |5T|) and 
recalling that (t| — Na^ 2Nap<^p^, and 2p(T^^^2 = —B we can now relate the upper left elements 
of both sides 

2Nap^p,-B = A{a,L;p,p) (52) 

justifying the earlier ad hoc decomposition of ^p^ (Sec. |8.2D and clarifying its basis. We see that 
variance sources {^pj and covariance sources (B) may both correspond to net two-point correlations 
(A). is a 'part' of A and not the reverse. This representation only hints at the richness 

of the system of scale-local correlation measures, even when restricted to pair correlations and 
(co)variance. Such structures also exist, albeit with increasing algebraic complexity, for arbitrary 
g-tuples and associated higher moments. 

12 Experimental Applications 

The previous material is quite general and is intended to accommodate at least conceptually the 
transition of any dynamical system across a phase boundary near which the dynamical DoF are 
ill-defined. However, experimentally we must contend with the hadronic final state of HI collisions 
and the measures that can be formed there. We now consider some practical measures and analyses, 
with examples recently proposed in the literature. 

12.1 Available measures and analysis types in the hadronic final state 

In the hadronic final state the basic event-wise measures rh are particle abundances and momenta, 
nij{6x) and pij{Sx), where i refers to particle species and j indicates a bin index. Special cases 
include nij{a) = 1, nij{L) = Ni ~ event abundance for species i, Pij{a) = particle momentum 
and Pij{L) = Pi - event total momentum for species i. These constitute the primary measures in 
fh = {mjj}. We can also add an integer index n randomly assigned to particles. 

A secondary measure system is derived from correlation analysis of individual events. Distri- 
bution statistics (< pt >) and model parameters (slope parameter) are simple examples. HBT 
correlation analysis in principle provides event-wise estimates of source radii, fiow analysis can 
estimate the reaction plane and azimuthal anisotropies, and more general event-wise correlation 
analysis can provide a complete topological assessment of each event in terms of information or 
dimension transport. A minimal set of correlation parameters which completely characterizes the 
event can then be included in fh. 

Depending on the abundance of a given particle species one can pursue an integral or differential 
analysis. In the case of limited abundance one can compare total or per-bin covariance matrix 
elements at two scales, element by element. These are just central-limit tests as represented by 
Eq. dll). In the case of greater abundance one can examine directly the integrands of elements of 
A((5xi, (5x2; m ?7^). One can analyze differences P2,obj{Sxi,6x2;mi,m2) — p2,refiSxi, 6x2;mi,m2) 
or preferably density ratios P2,obj / P2,ref to search for departures from a reference. 
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12.2 Partition systems 

The partition systems on P space directly accessible in the hadronic final state are at the particle 
(a) and event (L) scales. Additional partititions can only be defined on {m, n} and therefore only 
indirectly qualify as partitions on scale, to the extent that there are correlations between elements 
of {m, n} and P space. For instance, strong correlation between longitudinal momentum (rapidity) 
and axial position means that a partition on rapidity or pseudorapidity provides some ability to 
partition P space at different scales on axial position. However the scale interval is limited because 
this correlation is not perfect due to thermal and other fluctuations in longitudinal momentum 
about the basic Bjorken expansion correlation. Results obtained with such a partition must be 
carefully interpreted over a restricted scale interval. 

A partition on azimuth angle is also possible, but its interpretation depends on the degree of 
correlation between azimuth angle in momentum space and configuration space. Depending on the 
emission model this could be very stong (emission from an opaque source) or negligible (emission 
from a transparent source). Partition of the random index n is the most indirect but also least 
biased partition method, but its relationship to a scaled partition is not clear in general. 

12.3 Nonstatistical fluctuations 

The distinction between nonstatistical ('dynamical') and statistical fluctuations or variance can be 
seen as insubstantial given the discussion of scale-dependent variance presented above. 'Statistical' 
variance is accumulated below and within some scale interval characteristic of a correlation onset, 
or below the scale interval of interest. Dynamical fluctuations as discussed here simply correspond 
to variance accumulated on a complementary (larger) scale interval - the one of immediate interest. 
In a different analysis some statistical fluctuations might be considered dynamical by the same ar- 
gument. For discussion purposes we adopt a distinction of convenience between statistical variance 
and its complement on scale. 

The term A in Eq. (p2|), which represents the most global aspect of nonstatistical fluctuations, 
can be decomposed into several contributions as follows 

Aia,L;p,p) « {N {q,,, - 1 + aU + N{N - 1) {q^, - 1 + a^J} 

+ iV(iV-l)4 + ... (53) 

where ijres £ [1;2-|-] represents the effective sibling number for secondary particles derived from a 
single primary parent {e.g., resonance decays), qqc G [0,2] represents the effective sibling number 
for two-point quantum correlations (HBT and FSI), represents event-wise fluctuations in 

Qres, <7|^^ relates similiarly to qqc and o"?^ represents event-wise variation of the parent mean {e.g., 
distortion of the parent pt distribution by event-wise fluctuating DCC or radial flow contributions). 
Thus, possible EbyE fluctuations in single-point and two-point aspects of the parent distribution 
are here included with possible sources of two-point correlation as contributions to A. Mean values 
of quadratic forms are essentially the only aspects of EbyE variation of the parent distrbution 
accessible to a two-point or variance analysis. Note that an uncorrelated system {q = 1) and 
stationary parent (u^ = 0) result in j4 ^ 0, consistent with the CLT. This gives some impression 
of the complex task of interpreting a nonzero or zero value for A. Some unfolding of these possibly 
bipolar contributions can be achieved by noting multiplicity dependence, but there is not a unique 
inversion in general. 
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12.4 The subevent method 



In the subevent method [|T9[ each event is partitioned into two parts at large scale (~ L/2), for 
instance by a pair of disjoint bins on rapidity, azimuth angle or a random particle index n. One 
then extracts the covariance of one or a pair of additive measures m (= n{L/2),p{L/2) ■ ■ ■) or 
ratios <m>= mi/m2 obtained from the two parts. Effectively, the two-point density analyzed is 
p{L/2, L;mi,m2) on G2. We can interpret this method within the present context by referring to 
Eq. (|36|). The LHS is the variance difference referred to in |jl^ as the 'direct' method (in contrast 
to the subevent method). For specific cases (scale intervals) the RHS is the same covariance which 
forms the basis for the subevent method. For these cases the quantity A{5xi,6x2]Xi,X2) can be 



identified with the quantity dyn [|19i ) thus also providing the relationship between and the 
subevent method. (The discussion of in assumes lack of correlation between <pt> and N, 
so the corresponding covariance term B in Eq. (^) is missing.) 

One must carefully specify the scale interval being used for a given analysis. The subevent 
method as defined refers approximately to the scale interval [L/2,L], whereas the 'direct' method 
as discussed in [^] and the ^p^ measure as it is usually presented refer to the inteval [a,L]. The 
subevent method should therfore nominally be sensitive only to correlation developed or total vari- 
ance accumulated in the interval [L/2,L], and should not be directly comparable with an analysis 
based on [a,L]. However, because the partition methods employed are not applied directly to the 
primary space but only indirectly to a measure space (momentum space or a random index) it is 
possible for small-scale correlations to be included at some unknown level, prompting a recommen- 
dation to partition momentum space at large scale in such a way as to exclude issues of two-track 
resolution, BE/FD correlations and FSI. One effectively cuts on density p{a, L;pi,p2) in 6*2(^1,^2) 
so as to exclude a neighborhood around the main diagonal (at small scale) which contains quantum 
correlations and final-state interactions. 

Various multiplicity dependencies are considered in |19| as means to unravel multiple contribu- 
tions to 

dyn similar to the contributions included in Eq. ( |53| ) . Missing from |19|] are possible 
contributions represented by terms ct| for resonances and quantum correlations from event-wise 
variation of the two-point parent distribution. As with all global-variables analyses direct analysis 
of two-point density ratios on space G2(mi,m2) provides a more complete way to unfold differing 
contributions to variance comparison measures when the data have sufficient statistical power to 
support this approach. 

The subevent method is analogous to forming a mixed-pair reference from similar events, execept 
that the 'events' compared are subevent components of a single event. Use of large-scale event 
partitions to establish a statistical reference may be difficult to interpret in some cases, albeit 
computationally convenient. Assumptions must be made about the correlation structure of events 
which may be difficult to justify. An alternative for the study of large-scale correlations may be 
to compare sibling and mixed subevents - subevent pairs formed in similar ways from the same 
event and from similar events. This procedure would form part of a general scale-local approach 
to two-point correlations spanning the scale interval from particle pairs to subevent pairs, with the 
same general analysis method used throughout. 

12.5 Ratio fluctuations and resonance production 



In 1 30 1 it is proposed to study event- wise fluctuations in hadron multiplicity ratios in order to 
determine resonance abundances at chemical freezeout. The method is based on a change in the 
fluctuation structure of +/— multiplicity ratios to the extent that resonance decays (R) as well 
as thermal production (T) contribute to hadron yields. If a +/— multiplicity ratio is defined by 
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Ri2 = N1/N2 and we assume no correlations among resonances {a^ ~ A^,.) then the relative variance 
in the ratio is given by (r is a resonance type index) 

where 

91 = ^{4,1 + E^'-^} 

92 = ^{4,2 + 

512 = } Wt,12 +'^^r nr,l 71^,2} (55) 

vAlA'2 r 

If we further assume that Ni ^ N2 ^ N = Nt + Nr, Nr = J2r ^r, = Nr/N, ar = 1 - 
aji, (T^^j Nt, 4 12 ^ 0' = 1' then we have 51 ~ 52 ~ 1 and 512 ~ ctR. Therefore, 
^Ri2obj/'^Ri2ref ~ 1 " "iJ- Thesc assumptions apply to correlation-free thermal pion emission 
combined with a single thermal resonance type decaying strictly to tt"*" — 7r~ pairs as a simple 
example. 

Eq. ( |5^ ) for Gi(ni,n2) is a variant of Eq. (p3[ ) for Gi{p,n) (equivalent in form if the latter 
is divided through by P). There are two similar cases described by this algebra: distributions on 
(ni/n2,ni • 77-2) and on (ni — 712, ni + 712). The relative fluctuation distribution is essentially the 
same in the two spaces if a gaussian model and linearization are appropriate. The obj/ref ratio 
for multiplicity products (or sums) in the same case would be cj^jj obj/^ni2 ref ~ ^ + '^R [|l6|, where 
ni2 = Ni ■ N2. For the simple example considered, the thermal contribution is uncorrelated on 
Gi(ni,n2), with equal variance on all variable combinations, whereas the resonance contribution is 
perfectly correlated along Ni = N2. A superposition of these two contributions with fractions qt 
and aR gives the stated result. 

While correlations in the resonance contribution decrease the variance of the difference and 
ratio distributions they increase the variance of the marginal distributions on Ni and A2 and the 
sum and product distributions. The variance of < pt >= Pt/N is also increased by the presence 
of resoncances because the effecive sample number N^jf is reduced by correlations. In the limit 
of all resonances decaying to pion pairs and no thermal pions the observed number of particles is 
twice the number of statistical DoF as reflected in an effective sibling index q^es ~ 2 in Eq. (^3|), 
leading to a doubling of observed <pt> variance over the CLT expectation for the observed particle 
multiplicity. 

By some combination of fluctuation measurements of ratios/products or sums/differences and 
invocation of known resonance-decay schemes one can hope to infer the fractional resonance yield at 
chemical freezeout. Precision studies of resonance production, which are interesting in themselves 
and essential as a baseline to understand other aspects of the hadronic final state, require the 
formation of a reference system, so that one is dealing not with the ratio/product or sum/difference 
covariance system for real events alone but rather with differential departures from a reference. 
There is also the problem that other correlation sources may compete with resonance decay to 
alter A_|-/A_ fluctuations. 
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12.6 Covariance analysis of global variables 



In |31] it is proposed to increase sensitivity to possible QGP 'signals' by a combined analysis of 
the covariance of two or more global event properties. Event- wise antiproton multiplicity and 
transverse emitter radius are considered as an example. This approach is similar to the analysis 



format presented in Sec. |11.3| , but without use of a reference to form a comparison measure. 
The centrality dependence of covariance is itself taken as the indicator of unexpected behavior. 
One interesting aspect of this proposed application is that one of the measures is not a simple 
disttribution statistic like total multiplicity or mean pt but a more complex correlation aspect 
of the two-point momentum disttribution. This proposal is consistent with a general program 
to measure fluctuations of all kinds in the neighborhood of a phase boundary in order to better 
understand the underlying Lagrangian by enumerating and identifying the relevant DoF for various 
constraints. 



The underlying model premise in |31] is that, depending on centrality, collisions may fluctuate 
between two disconnected event classes, one {q) manifesting characteristics associated with QGP 
formation {e.g., increased emitter radius and antiproton abundance mean values) and the other 
[h) manifesting conventional hadronic systematic variation with participant number. The measure 
system in this case is (mi, 7712) = {Np,Rt), and the centrality dependence of the covariance is 
contrasted with that of the mean values of the measures for sensitivity in terms of this fluctuation 
model. While the paper emphasizes covariance measures it admits that variances would exhibit 
similar centrality dependence for this model. The central point is that attention to fluctuation 
measures is critical for detection of intermittent phenomena in collisions. 

The basic elements of the simulation are the estimated baseline measure fluctuations in the 
hadronic and QGP regimes, and the assumed shifts in mean values from one event class to the 
other. These assumptions determine all relevant aspects of the simulation. Mean values for the 
two event classes [h,q] are given by Np{b) [56, 78] {Nch{b)/Nmax) (antiproton multiplicity into 
STAR acceptance) and Rt{b) [6,9 fm] {Nch{b)/Nmax)^^'^ , where N^h is the charged-particle 
total multiplicity used to infer event centrality and N^ax is the mean value for central collisions. 
Variances are approximated by a'j^_ « Np and o"^^ « B^, the second dominated by an estimate of 
measurement error rather than intrinsic statistical fluctuations. The result is a rather large increase 
in covariance in a crossover region of centrality where the two event types are nearly equiprobable. 
We observe in a scatter plot on {Np,Rt) that the joint distribution is bimodal due to the chosen 
relationship between assumed mean-value shifts and variances. 

The assumption of catastrophic fluctuations between two event types maximizes the effect of 
the model on variances and covariance. An alterative scenario consisting of a partial transition (still 
first order) over some fraction of the event volume (having its own fluctuation structure) would 
produce a significantly smaller effect on correlation measures. These model assumptions do however 
correspond to some popular QGP scenarios. We note that an expected major source of covariance, 
collision geometry and stopping fluctuations ('volume' fluctuations) as illustrated in Fig. |^ and 
present for either event type, is here probabably underestimated by the Poisson variance of N^h- 

Having modelled the most dramatic QGP scenarios one might then consider how to establish 
by these correlation measures either that a QGP hypothesis is falsified at some level, or that some 
fraction of some events does indeed undergo a transition to a nonhadronic phase. In the adopted 
scenario the input to the corelation measures is bimodal, a rather strong statistical manifestation. 
A weaker and possibly more realistic situation to simulate is one in which anomalous events form a 
tail structure on a hadronic core. In this case one has to compare carefully the centrality dependence 
of covariance matrix elements to that for a reference and perform a significance study of the result. 
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There is much technical work needed to elaborate these and other EbyE analysis methods. We 
see from these examples and interpretations that a general description of scale-local covariances 
and comparison measures permits a comprehensive interpretation of some apparently dissimilar 
methods. 

13 Toward a Generalized Central Limit Theorem 

From the original concept of the elementary central limit theorem as an integral variance compar- 
ison measure for a gaussian fluctuation model, we have enlarged the analysis to multiple measures 
described by a covariance matrix, have expressed the basic comparison in a scale-dependent man- 
ner, have espressed the source of departure from the CLT in terms of a matrix of scale integrals 
comparing object and reference two-point densities and have considered as the most differential 
form of comparison a direct analysis of density ratios on two-point measure spaces. 

From the vantage point of this scale-local picture we can now generalize the central limit theorem 
and reinterpret its role in the description of many-body systems. In terms of the new concepts the 
elementary CLT as presented in Sec. ^ compares total variance for a single measure at two points 
on scale to determine whether additional variance has accumulated over the interval. The CLT 
states that if there is no net correlation introduced within a scale interval the total variance of a 
measure is the same at each endpoint - T,'^{p,p, 8x2) — = 0. 

We have subsequently developed a greatly extended version of this elementary comparison 
procedure and generalized it to a pair-wise total (co)variance comparison for an arbitrary measure 
system as represented by 



and have given the general form of the source of any total variance change in terms of a matrix of 
covariance integrals over two-point measure densities represented by A. We can therefore generalize 
the CLT to a statement about the scale invariance of covariance: 

The total covariance matrix IC{m,5x) for a measure system on the space Q is scale 
invariant over a bounded scale interval [Sxi,5x2] in the absence of all categories of 
net two-point measure correlation over the same scale interval. Departure from scale 
invariance is represented by covariance integrals of density differences on corresponding 
two-point measure spaces. 

In this formulation we have represented the possible sources of accumulated covariance as integrals 
over two-point measure spaces. We have also introduced the possibility that a uniform reference 
(implicit in standard variance calculations) may be replaced by an arbitrary reference or model 
appropriate to a given problem. The resulting system offers the opportunity to choose from a 
continuum of techniques and references, depending on the fluctuation model and statistical power 
of the data. 

14 Conclusions 

We have as our nominal goal identification and study of a quark-gluon plasma, a simplistic ter- 
minology for a complex phenomenon. It may be that the emergence of a color-deconfined QCD 
system is subtle and topological rather than catastrophic and thermodynamic. Production of a 
QCD plasma as a catastrophic process with discontinuous signals may not be forthcoming. To 




(56) 
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insure the best opportunity for discovery we require sensitive quantitative comparison measures 
which extract complete correlation information from collision events. 

In studying the correlation structure of HI collision events at the CERN SPS we have been struck 
by their sameness and ordinariness, the apparent close adherence to simple statistical fluctuations 
in the distribution of mean pt and the wide-ranging description of hadronic ratios and abundances 
by a simple statistical model with a universal chemical freezeout temperature. These two results 
are deeply connected by their relationship to the central limit theorem and related maximum 
symmetry principle. A common supposition has been promoted that QCD plasma formation at high 
energy densities may lead to dramatic, or easily quantifiable, changes in the correlation structure 
of the multiparticle final state, changes detectable by elementary correlation measures. Given the 
unspectacular SPS results we are now forced to re-evaluate our need for sensitivity. 

A central goal in EbyE analysis has been to develop global comparison measures sensitive to 
excess variance which might signal residual correlations due to incomplete equilibration - possibly 
structure remaining from a phase transition. We have seen that these global measures are ultimately 
based on the central limit theorem as a statistical reference, and the CLT is in turn based on a 
scale-integral concept. The present treatment makes it clear that global- variables analysis based 
on scale integrals - an unavoidable strategy in case of low statistical power - is rather limited in 
its sensitivity, interpretability and power to discriminate among various correlation sources. 

Statistical power (event multiplicity) determines the optimum scale interval and scale resolu- 
tion of a correlation analysis. Low statistical power requires the use of integrals over space and 
scale to extract statistically meaningful correlation measures, with accompanying loss of informa- 
tion. Higher statistical power makes possible extension of the analysis to scale-local differential 
measures and direct analysis of densities on two-point measure spaces. One can then study the 
scale dependence of density ratios to extract full correlation information from data. Future de- 
velopments include direct determination of two-point densities, analysis of scale-local topological 
measures extracted from these densities and extention of the analysis to higher moments {q> 2). 

The question of scale invariance represented by the CLT is deeply related to the nature of 
distribution symmetries. Scale-invariance is equivalent to constant autocorrelation - Eq. ( |38[ ) - 
which is in turn an indication of maximum symmetry. Thus, the CLT in its generalized form is a 
statement about symmetry variation over a bounded scale interval, a potentially new insight. With 
a uniform distribution as reference one can inquire whether an object distribution is maximally 
symmetric over a scale interval and how its symmetries evolve dynamically. More generally, one 
can ask whether object and reference distributions have the same symmetry or autocorrelation over 
some scale interval. This presents a new approach to quantitative description of nonequilibrium 
systems and the equilibration process. 

Scale-dependence of symmetry variation is the underlying issue in the study of phase transitions. 
Near the QCD phase boundary hadrons and partons as dynamical objects are transitory and ill 
defined. Basing a dynamical description on a particle hypothesis near a phase boundary is a poor 
strategy. To provide a fundamental treatment one must ultimately adopt a descriptive system in 
which particles are not fundamental elements, in which arbitrarily correlated measure systems play 
the central role. A generalized central limit theorem provides a basis for such an approach. One 
can study the symmetry dynamics of a measure system relative to a maximally symmetric system 
or other symmetry model. 

In a broader sense we are probing the dynamics of vacuum symmetries within a violent collision. 
The more traditional study of symmetry statics is carried out in asymptotic regions of a constraint 
space far from any phase boundary. Symmetry dynamics is important in the neighborhood of a 
phase boundary. The QCD phase boundary marks a region of constraint space where the nature of 
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vacuum symmetries is changing dramatically and fundamentally. This is a challenging undertaking 
requiring the best possible analysis techniques. 
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